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Appendix  B 


Section  1 


INTRODUCTION  AND  SUMMARY 


This  is  the  final  report  on  AFGL  Contract  F19628-81-C-0048  titled 
"Electron  Beam  Transport  in  the  Ionosphere  -  Energy  Deposition  and  Optical 
Emissions  Based  Upon  the  Combined  Effects  of  Plasma  Turbulence  and  Particle- 
Particle  Interactions".  This  work  has  been  carried  out  for  the  Air  Force 
in  support  of  planned  rocket  and  satellite  kilovolt  electron  beam  (e-beam) 
experiments.  Our  goals  in  this  program  have  been  to 

1)  model  plasma  turbulence  as  it  affects  kilovolt  e-beams  in 
the  E-region, 

and 

2)  make  quantitative  estimates  of  beam  energy  deposition  and 
luminosity  as  functions  of  altitude  z  and  radial  distance  r. 

This  work  is  an  extension  of  our  efforts  under  a  previous  contract 
in  which  code  MCBE  was  developed  and  applied  to  the  problem  of  particle-particle 
interactions.  MCBE  is  a  2D  Monte  Carlo  code  explicitly  treatinq  electron  motion 
in  a  magnetic  field.  Energy  loss  is  approximated  as  continuous  and  beam  spread¬ 
ing  comes  about  by  use  of  a  multiple-scattering  formula. 

Particle-particl e  interactions  are  not  sufficient  to  cause  the 
spreading  observed  in  some  beam  produced  luminosity  streaks  of  past  experiments. 
(Davis  et  al.  (1971)  and  Hallinan  et  al.  (1978)).  We  were  able  to  determine  the 
extent  of  this  deficiency  in  our  earlier  work  using  MCBE.  It  is  clear  from  this 
work  as  well  as  from  an  examination  of  Hallinan's  data  that  plasma  turbulence  can 
play  an  important  role  in  e-beam  propagation  in  the  E-region.  Much  of  our  effort 
under  the  current  contract  has  been  to  quantify  the  effect  of  this  turbulence. 
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It  is  helpful  to  briefly  outline  the  approach  we  have  taken  in 
modeling  beam-plasma  interactions  prior  to  presenting  the  details  in  Sections 
2  and  3.  We  first  assume  the  wave-particle  and  particle-particle  regions  can 
be  effectively  separated.  The  basis  of  the  model  is  then  that  the  given  tur¬ 
bulence  can  cause  a  series  of  effective  displacements  of  the  beam  electron  guid¬ 
ing  centers  through  a  corresponding  series  of  briefly  experienced  perpendicular 
components  of  the  turbulent  electric  field,  Ex  .  The  lower  altitude  boundary 
for  this  effect  is  defined  as  the  altitude  at  which  the  turbulent  growth  rate 
egua I s  the  particle-particle  collision  frequency.  Above  this  altitude,  we  use 
a  Monte  Carlo  description  to  follow  the  motions  of  the  beam  electron  guiding 
centers  caused  by  the  above  mentioned  displacements  if  they  are  given  altitude 
dependence.  Otherwise,  the  final  result  of  the  displacements  is  given  by  an 
analytical  radial  distribution  function  (Gaussian). 

A  given  displacement  is  a  function  of  both  Ex  and  an  auto¬ 
correlation  time  ;.  We  choose  to  replace  Ex  by  an  expression  containing  ;v, 
the  velocity  spread  of  the  beam  reflecting  the  strength  of  the  turbulence. 

The  nature  of  this  quantity  will  be  considered  for  both  cold  and  warm  beams. 
Emphasis  will  be  on  warm  beams  for  which  sv  can  be  equated  to  the  beam  thermal 
velocity.  Specification  of  6v  for  this  case  is  made  through  the  parameter 

O 

.  which  gives  the  ratio  of  the  thermal  to  the  beam  energy  [,:  =  (-v/v)  ]. 
Spreading  is  then  examined  for  ranging  from  .01  to  .1  with  the  guiding 
center  displacement  being  constant  versus  altitude.  It  is  difficult  to  pre¬ 
scribe  the  beam  thermal  velocity.  For  this  reason  we  have  chosen  limiting 
values  of  ;■  for  a  given  beam  so  as  to  bound  the  beam  spread.  This  scoping 
approach  to  the  problem  is  further  justified  because  of  the  large  uncertainties 
in  the  experimental  measurements. 

Following  the  above  introductory  remarks,  we  summarize  our  results 

as  follows: 

1)  A  model  has  been  developed  which  provides  for  beam  spreading 
by  plasma  turbulence.  The  model  is  based  on  random  walk 
motion  of  the  beam  electron  guiding  centers  caused  by  the 
perpendicular  component  of  the  turbulent  electric  field. 


2)  The  model  has  been  used  in  an  analysis  of  beam  spreading 
reported  by  Hallinan  et  al.  (1978)  [ECHO  IV  experiment]  and 
Davis  et  al.  (1971).  We  can  account  for  most  beam  widths 
by  adding  the  additional  spreading  due  to  particle-particle 
effects.  For  the  pulses  injected  below  180  km,  our  fits 
suggest  a  [,  value  <  .01. 

3)  Detailed  energy  deposition  results  (in  r  and  z)  are  pre¬ 
sented  for  two  of  the  beam  cases  considered.  The  results 
come  from  code  MCBE  for  incident  radial  distributions  given 
by  a  delta  function  and  by  the  above  mentioned  turbulence 
model . 

4)  Optical  properties  of  the  beams  treated  in  3)  are  determined. 
These  include  r,z  distributions  of  surface  brightness  for  the 
sum  of  most  emission  features  in  the  LIV  and  visible  portions 
of  the  optical  spectrum.  Scaling  factors  are  also  provided 
which  relate  emission  efficiencies  to  energy  deposition  for 
the  individual  features  just  mentioned. 

As  a  concluding  remark,  the  work  carried  out  in  this  program 
together  with  that  from  the  earlier  one  enables  us  to  examine  the  r,z  dependence 
of  energy  deposition  and  optical  emissions  for  electron  beams  propagating 
through  the  E-region  under  the  influence  of  particle-particle  interactions 
and  plasma  turbulence. 
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Section  ?. 


RANDOM  WALK  MODEL  FOR  WAVE-PARTICLE  INTERACTIONS 


In  this  section,  we  will  describe  the  random  walk  of  the  guid¬ 
ing  center  of  the  beam  electrons  due  to  wave-particle  interactions.  This 
type  of  treatment  is  meant  to  place  the  mechanism  of  the  spreading  caused 
by  wave-particl e  interactions  on  the  same  footing  as  that  for  particle- 
particle  interactions;  the  details  of  the  latter  have  been  given  by  Lin 
and  Strickland  (1981).  In  the  case  of  particle-particle  interactions, 
the  random  force  is  provided  by  collisions  between  beam  electrons  and 
neutral  particles.  For  wave-particle  interactions,  the  random  force  is 
due  to  the  turbulent  E  field  acting  on  the  beam  electrons. 


2 . 1  Guiding  Center  Motion . 

Let  Ef  be  the  magnetic  field  ana  z  the  direction  of  5. 
The  equation  of  motion  of  a  beam  electron  is  given  by 


dy  _  ev  x  §  ef 
dt  ~  me  m 


We  are  not  interested  in  the  coord;nates  of  the  beam  electron; 
instead,  what  concerns  us  are  the  coordinates  of  its  guiding  center.  In  the 
x-y  plane,  the  electron  is  positioned  at  (x,y)  with  velocity  (v  ,v  )  while 

x  y 

its  guiding  center  is  at  r^  =  ( XGC ’^GC ^ aro  rp^ated  hy 


xrr  =  x  -  v  / 
GC  y 


yGC  =  Y  +  V 


(2) 
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where  '2  is  the  gyrof requency. 


Taking  the  perpendicular  component  (i.e.,  the  component  in 
the  x-y  plane)  of  Equation  (1)  and  using  Equations  (2)  and  (3),  we  arrive 
at 


or 


fee  = 

d  t  m 


Ar 


GC 


(4) 


(5) 


where  E^  is  now  the  perpendicular  component  of  the  turbulent  E  field  and 
r  is  the  time  interval  for  a  coherent  interaction  (the  auto-correl ation 
time) . 


Equation  (5)  is  the  basis  for  our  random  walk  model  for  wave 
particle  interactions.  It  can  also  be  derived  using  the  E*X[f  drift  velocity 
?X§/cB2  to  get  Ar  =  v^^  t.  The  picture  of  the  random  walk  of  the  guiding 
center  is  described  in  the  following  subsection. 

2 . 2  Spatial  Properties  of  the  Turbulence  Model. 

The  turbulent  plasma  is  considered  to  contain  regions  of 
appreciable  electric  field  strengths.  Within  each  region,  phase  coherence 
is  maintained  in  the  manner  described  by  Stix  (1964).  The  phase  of  oscil¬ 
lations  in  one  region  is  random  with  respect  to  the  phase  in  each  other 
region  due  to  nonlinear  effects  of  the  large-amplitude  disturbance.  The 
length  within  which  coherence  is  maintained  is  called  the  correlation  length. 
(Stix  ' s  analysis  results  in  a  correlation  length  of  about  0.64  '  where 
is  the  effective  wavelength  of  the  turbulence). 

In  a  correlation  length,  the  beam  electron  sees  a  constant  E 
field  (see  Appendix  B  for  details).  Each  of  the  three  dimensions  of  a 
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coherent  region  is  about  one  correlation  length.  Thus,  the  beam  electrons 
feel  the  constant  E  field  during  the  time  they  traverse  a  correlation  length 
and  then  run  into  another  coherent  region  where  they  feel  another  constant 
t-field  whose  phase  is  random  with  respect  to  the  E  field  in  the  previous 
coherent  region,  and  so  on.  Ine  guiding  center  accordingly  undergoes  a  random 
walk  motion  such  that  the  beam  spreads. 

2.3  Displacement  Ar^  as  a  function  of  Turbulent  Velocity  !v. 

We  first  define  needed  quantities  associated  with  the  coherency  of  the 
turbulent  spectrum.  The  correlation  length  2.  is  given  by 


=  ciA 


(6) 


and  the  coherency  time  by 


i  =  2/vb  (7) 

The  above  expressions  are  appropriate  for  coherent  wave  packets  with  group 
velocities  much  smaller  than  the  beam  velocity.  The  value  of  is  of  order 
unity  with  its  specific  value  given  by  the  applied  turbulence  model. 

Tne  potential  energy  of  the  turbulent  electrostatic  mode  can  be 

wri tten 


e*  =  8  m(<Sv)2]  (8) 

where  the  velocity  v  is  an  equivalent  measure  of  the  strenqth  of  the 
turbulence.  The  numerical  factor  in  Equation  (8)  corresoonds  to  the  re¬ 
sult  of  Drummond  et  al.  (1970).  For  cold  beam-weak  turbulence,  v  =  v,  -  /k 

u||  II 

as  discussed  by  Drummond  et  al .  (1970)  and  O'Neil  et  al  (1971).  For  strong 
turbulence,  ;v  •-  v^/4  ,2  as  described  by  Stix  (1964).  We  assume  the  cold  beam-weak 
turbulence  form  can  also  be  applied  to  the  warm  beam-weak  turbulence  case  which  shall 
be  of  particular  interest  to  us.  In  Section  3.1,  this  will  lead  to  v  being 
equated  to  the  beam  thermal  velocity. 
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The  wave-particle  interaction  is  a  collisionless  effect.  Mode 
damping  will  occur  when  particle-particle  interactions  become  large  enough 
at  the  lower  altitudes.  This  limiting  particl e-particle  frequency  can  be 
wri tten 

(particle-particle  collision)  =  C„,  (9) 

where  r  is  the  linear  growth  and  is  of  order  unity. 

Let  be  the  angle  between  k  and  £  and  b  the  angle  between 
vb  and  5.  We  introduce  an  effective  gyroradius  r  : 


r 


g 


*  V 


(10) 


which  is  defined  independent  of  the  angle  . 

With  the  use  of  Equations  (6),  (7),  (8),  and  (10)  we  can 
write  Equation  (5)  in  the  form 

01^-  *  8  /FrCj  sin  (fy  (11) 

9  b 

Use  is  made  of  E_,_  =  E|(  tan  3^,  qk)(  =  E(|  ,  and  kA  =  2'.  To  apply 
Equation  (11)  to  a  specific  problem,  we  must  provide  a  specification  to 
v  which  will  be  done  in  Section  3. 


2.4  Ra_dJ_a_l_  Distribution  Function  and  Avearge  Beam  Radius. 


In  this  subsection,  we  discuss  the  two  cases  where 


'  GC 


1  s 


constant  and  variable.  For  the  case  of  constant  Ar^ ,  the  needed  function 
for  the  beam  radial  distribution  after  N  displacements  of  the  guiding 
center  is  analytical.  For  the  other  case,  a  numerical  procedure  is  called 
for-  and  this  has  been  developed  by  us  during  this  program.  Dotailr  now  follow. 
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We  begin  by  defining  N,  the  total  number  of  steps  or  guiding 
center  displacements  caused  by  the  turbulence.  It  is 


N 


{12) 


where  t  is  the  travel  time  of  the  beam  electrons  within  the  region  where 
the  wavc-particl e  interaction  is  dominant.  Since  the  parallel  wave  pnase 
velocity  is  equal  to  the  corresponding  component  of  the  beam  velocity  in  the 
types  of  plasma  turbulence  to  be  considered,  viz. 


COS  t 


'b  ’ 


(13) 


we  can  write  the  coherence  time  in  terms  of  the  wave  frequency: 


=  2r<  Cj  cos  cos 


(14) 


If  we  let  d  be  the  distance  between  the  injection  altitude  and  the  altitude 
at  which  the  particl e-pa rticl e  collision  frequency  is  C?  times  the  linear  growth 
rate,  then  t  =  d/V^  cos  ^  and  we  obtain 


N 


_ 1 

Cl2~  cos  k 


I 


IS) 


A  typical  case  of  =  6  x  10  sec  ,  d  =  50  km,  mv,/2  =  26  keV,  L  =  45  , 
0  b  k 

and  ^  -  70u  generates  an  N  of  8512  if  Cj  =  1. 


H 


Let  W^(r)  designate  the  radial  distribution  function  after  N 
displacements  where  N  refers  to  an  average  for  each  of  the  beam  particles. 
For  constant  Ar^,  W^(r)  has  the  form 


with  normalization 


WN(r)  2i;rdr  =  1  .  (17) 

W^(r)  is  the  same  distribution  function  given  by  Lin  and  Strickland  (1981)  for 
particle-particle  interactions,  except  now  Ar^  is  to  be  specified  for  wave- 
particle  effects. 


We  may  go  one  step  further  for  the  constant  step  size  case  and 
express  analytically  the  average  radius  of  the  beam  from  Equation  (16). 
Designating  it  by  <r>w  (in  units  of  rg),  it  is 


<r 

w 


J<  ./^T  Av_\2 
C0S  -b  V  vb  ) 


(IE) 


We  will  make  use  of  tquations  (16)  and  (18)  in  the  next  section 
as  we  discuss  the  ECHO  IV  experiment. 
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For  non-constant  A r^ ,  we  have  developed  a  Monte  Carlo  code(WPDIF) 
to  describe  the  2D  random  walking  of  the  beam  electron  guiding  centers. 

A  brief  description  of  the  input  parameters  and  calculation  follows  starting 
with  the  parameters.  They  are: 

Zq  starting  al ti tude 

eg  cosine  of  the  pitch  angle  of  the  starting  beam  electrons 
E  electron  energy 

n^(z)  altitude  profile  of  the  plasma  density 

I.  starting  beam  current 

fg  starting  distribution  of  guiding  centers 

M  number  of  histories  or  electrons  to  be  sampled. 

The  value  of  M  used  to  obtain  our  initial  results  is  1000.  These  M 
electrons  are  allowed  to  advance  in  their  random  motion  as  a  group.  Thus 
in  passing  through  a  given  A z  corresponding  to  the  coherence  time  ~ ,  M 
guiding  center  displacements  are  determined  and  recorded.  This  enables 
us  to  specify  the  decreasing  beam  density  at  each  step  for  input  into  the 
expression  for  the  displacement  through  the  next  step.  The  calculation 
proceeds  until  the  growth  rate  y  for  the  given  instability  is  exceeded  by 
the  electron-neutral  collision  frequency. 

This  completes  the  formalism  of  our  random  walk  model.  We 
now  turn  to  applications  in  which  specification  will  be  given  to  ov  and  N. 
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APPLICATIONS  OF  THE  RANDOM  WALK  MODEL 


We  now  consider  the  level  of  plasma  turbulence  anticipated  for 
beam  parameters  of  interest.  It  is  convenient  to  follow  Davidson  (1972)  in 
defining  two  categories  of  beam  driven  turbulence.  Weak  turbulence  corresponds 
to  turbulent  energy  densities  much  smaller  than  the  kinetic  energy  of  the  beam, 
whereas  strong  turbulence  implies  that  a  large  portion  of  the  beam  energy  is 
transferred  to  turbulence.  Since  beam  energy  losses  to  plasma  turbulence  have 
been  found  experimentally  to  be  small,  we  consider  weak  turbulence  in  this 
report . 


We  distinguish  two  classes  of  weak  turbulence  by  examining  the 

linear  theory  of  beam  plasma  instability.  The  growth  rate  depends  on  the  beam 

temperature;  if  the  temperature  is  sufficiently  large,  the  condition  for  a 

1/3 

warm  beam  is  V|3tf/Vb  >  a  ,  where  vbtb  is  the  beam  thermal  velocity  and  >  is 
the  ratio  of  the  beam  to  ambient  electron  density.  When  the  inequality  is  re¬ 
versed,  the  beam  is  considered  cold.  It  is  difficult  to  access  the  values  of 
the  beam  temperature.  For  cases  where  a  is  quite  small,  however,  we  assume 
the  beam  temperature  can  satisfy  the  above  warm  beam  criterion.  In  what  fol¬ 
lows  we  apply  the  results  of  warm  beam  theory  to  the  high  energy- low  current,  and, 
thus,  small  a  beams  characteri sti c  of  the  ECHO  IV  experiment  (Hallinan  et  al . 
(1978)).  On  the  other  hand,  we  use  cold  beam  theory  to  explain  beam  spreading 
in  the  rocket  experiment  reported  by  Davis  et  al.  (1971)  where  a  relatively  low 
energy  (8.7  keV)  high  current  (500  mA)  beam  was  injected. 

3.1  Warm  e-beam. 


For  warm  beams,  >-v ,  is  a  constant  and  is  equal  to  the  beam  thermal 
velocity  vbtb.  This  may  be  seen  using  the  linearized  theory  of  the  two-stream 
instability  (see,  e.g.,  Stix  (1962)).  The  beam  and  the  background  plasma  con¬ 
sist  of  a  doubly-humped  magnetized  plasma  which  is  vulnerable  to  the  two-stream 
instability.  The  dispersion  relation  and  the  frequency  are  obtained  by  standard 


a 


methods  (Stix  (1964)).  T  he  wave  number  is  obtained  by  maximizing  g * b ( v  |(  ) 
which  is  the  derivative  of  the  beam  electron  distribution  function  for  the  zero- 
order  velocity  component  along  B.  The  form  of  is 

t7  exp  [-"rfv'b//2lT«b]  (19) 

The  above  maximizing  leads  to 


w/kii  =  V 


-  v 


bth 


(20) 


where  vfath  =  / kT„b/m.  Equation  (20)  indicates  that  the  parallel  beam  velocity  is 
slightly  greater  than  the  parallel  phase  velocity  by  the  amount  of  the  beam  thermal 
velocity.  Recalling  now  the  discussion  following  Equation  (8),  we  see  that 
Equation  (20)  leads  us  to  the  desired  result,  namely 


v  ~  vbth 


(21) 


The  magnitude  of  the  constant  step  size  is  probably  accurate 
during  the  initial  stage  of  the  beam  spreading.  This  may  not  be  as  good 
an  estimate  at  the  final  stage  where  the  random  walk  suddenly  stops.  To 
make  the  transition  smooth  from  the  region  where  wave-particle  interactions 
dominate  to  the  region  where  particle-particle  interactions  dominate,  we 
assume  that  the  overall  spreading  of  the  beam  is  effectively  the  same  as 
the  case  where  the  effective  constant  step  size  is  half  the  maximum  value 
given  by  Equation  (11). 

The  beam  radius  as  a  function  of  the  traveled  vertical  distance 
d  can  be  identified  with  the  mean  beam  radius.  From  Equation  (18)  we  have 


r> 


=  Aj  \  d 


Vh,2 

vb 


where  the  final  value  of  d  is  obtained  by  the  criterion  that  the  linear 
qrowth  rate  is  the  same  as  the  collision  frequency  (we  assume  C,,  =  1). 


(22) 
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It  should  be  noted  that  Aj ,  as  may  be  obtained  from  Equation  (18)  is  here 
to  be  halved  based  on  the  comments  of  the  previous  paragraph.  For  a  warm 
beam  the  linear  growth  rate  is  linearly  proportional  to  a  (Stix  ( 1 9b4 ) 
and  Tsvtovich  (1977)). 


■>  *  S  ■'  ■■  Q  (">  <«> 

bth 

where  Q(w),  a  function  discussed  in  Appendix  A,  is  evaluated  at  its 
maximum.  For  parameters  of  interest,  QTOX  ^  0.4.  The  important 

quantities  of  concern  here 'are  a  and  8.  y  is  expressed  as 

y  =  A2  ■./?■  (24) 

where  A?  contains  the  terms  in  Equation  (23)  other  than  «.  and  8  with  a  Q 
value  of  .4.  If  we  denote  the  initial  value  of  y  by  yQ,  then 


and 


y 


o 


Vfi 


(25) 


Y  =  y0/(1  +  R)2  =  Y0/(l  +  A1  /d  ?)2 


(26) 


where  we  have  assumed  the  thermal  energy  of  the  beam  remains  the  same  and 
the  beam  initially  has  a  radius  of  one  equivalent  gyroradius  r  . 

Equation  (26)  gives  the  linear  growth  rate  as  a  function  of  d. 
Final  d,  called  d^,  is  obtained  by  equating  y  to  the  collision  frequency. 

The  values  of  the  parameters  for  ECHO  IV  experiments  are  assumed 

as  follows: 
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C1  =  0.64 


ng  =  105  cm’^ 


107  sec  ^ 


The  quantities  ...  and  0^  can  be  worked  out  following  the  method 
prescribed  by  Stix  (1964).  The  resulting  values  are  ..  =  6  x  10^  sec"^  and 
=  45^.  Table  1  lists  the  values  of  parameters  for  all  observed  streaks  in 
ECHO  IV  experiments.  The  values  of  <jg  are  given  for  a  beam  density  distribu¬ 
ted  uniformly  out  to  the  initial  gyroradius  r  The  last  column 

2 

is  the  minimum  value  of  (vbth/vb)  satisfying  the  warm  beam  condition.  We 
see  in  all  cases  that  our  lower  limit  of  .01  for  B  exceeds  this  required 
minimum  value. 


Table  1.  Parameter  Values  of  ECHO  4  Experiments 


Pulse  Number 

Energy  (keV) 

ao 

°b 

min(vbth/vb)2 

12 

38 

.00023 

60° 

-  0C38 

10 

26 

.00035 

70° 

.0050 

21 

34 

.00036 

70° 

.0051 

3 

36 

.00018 

70' 

.0032 

24 

33 

.00042 

O 

c 

.0056 

1 

26 

.00027 

70° 

.0042 

The  decreasing  of  the  linear  growth  rate  and  the  increasing 
of  the  collision  frequency  for  the  experimental  parameters  are  shown  in  Figs.  1-3. 
Growth  rates  are  shown  for  the  two  B  values  of  .01  and  .1.  Two  collision 
frequencies  are  also  shown  reflecting  the  uncertainty  in  the  scattering 
cross  section.  This  cross  section  is  given  by  the  screened  Rutherford 
formula 
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ALTITUDE  (km) 


<r>w  FREQUENCY  (sec'1) 

Figure  1.  Growth  Rate,  Particle-Particle  Collision  Frequency, 
and  Mean  Beam  Radius  For  Warm  Beam  Model  Applied  to 
Two  of  The  ECHO  IV  e-Beams  Labeled  as  Pulses  1  and  10. 
Only  <r>w  for  pulse  10  is  shown  since  this  part  of  the 
figure  is  to  just  illustrate  the  general  behavior  of  the 
spreading  with  propagation  distance. 


ALTITUDE  (km) 


1 


ALTITUDE  (km) 


Figure  3.  Similar  to  Figure  1  Except  for  Pulses  3,  21,  and  2 4 
Again,  only  one  example  of  <r-  is  included  as  was 
the  case  in  Figure  1. 


(27) 


,E|  -  l' * 

eids  v2p2  ;rr-r+  it 

where  Z  is  the  atomic  number  of  the  scatterer,  v  and  p  are  the  velocity  and 
momentum  of  the  incident  electron,  and  r,  is  a  screening  parameter.  This 
latter  quantity  may  be  expressed  in  the  form 


4.3  Z2/_3 

’■""T"  'r 

where  •  is  constant  and  is  dependent  on  the  screening  potential.  Here, 
we  consider  two  values  of  g  (.6  and  1.0)  chosen  to  reflect  the  uncertainty 
in  this  potential  for  our  application  (Berger  et  al.  (1970)). 

It  is  interesting  to  observe  from  these  plots  that  the  cutoff 
altitudes  between  the  wave-particle  dominant  region  and  the  particle-particle 
dominant  region  are  about  150  km  for  all  the  ECHO  IV  experiments  if  ;  =  .01. 
If  -  .1,  the  cutoff  altitude  is  slightly  below  the  injection  altitude. 


The  range  of  final  beam  radius  for  P  ranging  from  .01  to  .1  is 
presented  in  Table  2  for  each  experiment.  The  n  value  of  .6  was  used  to 
specify  the  cut-off  altitude.  The  spread  for  pulse  12  is  narrower  than  that 
of  the  others  emu  is  due  to  a  combination  of  its  energy  and  injection  pitch 
angle. 


lable  2.  Beam  Spreading  Due  to  xlave-l'art  ii  I  e  i  ntei  <n.;  • .  >r. 


Pulse  Number 

Energy  (keV) 

<!'■ 

W 

df  (km) 

12 

38 

4.4  -  4.6 

.524  -  56 

10 

26 

4.6  -  6.6 

.22  -  45 

21 

34 

3  -  5 

.1-28 

3 

36 

.7  -  2.8 

.006  -  10 

24 

33 

1.3  -  3.4 

.02  -  14 

’* 

■>6 

.33  -  2.4 

.001  -  6  j 

. 

-  .  -  .  .  -  .  . 

1 

(28) 


Ip  Figure  4,  we  show  the  beam  spreading  due  to  particle- 
particle  interactions  (the  solid  vertical  line)  and  the  total  spreading  due 
to  both  particle-particle  and  wave-particle  interactions  (the  shaded 
bars),  txcept  for  the  26  keV  case  injected  at  210  km,  good  agreement  is 
found  for  all  the  other  cases  between  experiments  and  calculations. 


3.2  Cold  _e-J) earn 

We  now  consider  the  cold  beam  model  for  the  case  of  the  8.7  keV, 

500  m/\  beam  in  the  Davis  et  al .  (1971)  experiment.  Linear  theory,  giving  the 
dispersion  relation  for  cold  beams  injected  into  a  plasma,  predicts  (O'Neil  and 
Malmberg  (1968))  the  following  maximum  growth  rate 


) 


/f 

2 


ivl/3 

2' 


(29) 


where  u>  is  the  plasma  frequency. 

The  electric  field  corresponding  to  the  nonlinear  development  of  the  beam- 
plasma  instability  was  given  by  Drummond  et  al.  (1970)  and  O'Neil  et  al. 
(1971) : 


lc  _  1  /'ti  1/3  ,1  n  2  , 
8;;  ~  2  ^  ^2  b  b  ) 


(30) 


which  can  be  derived  by  the  usual  particle  trapping  and  bounce  frequency 
arguments.  For  this  case,  <Sv  is  given  by 


.‘v 


1 

2 


(31) 
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Figure  4.  Mean  Beam  Radii  For  ECHO  IV  e -Beams.  The  Vertical 
Line  is  From  Lin  and  Strickland  (1981)  and  Gives 
the  Approximate  Particle-Particl e  Contribution. 

The  Points  and  Error  Bars  Come  From  Hallinan  et  al. 
(1978).  The  Shaded  Regions  Give  the  Model  Results 
for  i-;  Ranging  From  .01  to  .1.  The  Particle-Particle 
Contribution  is  Included  in  These  Results. 
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Substituting  this  equation  into  Equation  (11),  we  get 


(,\r) 


GC  _  ,,  o* .  r  ■  ,  \2/ -i 

2  C |  sin  ^  ^ 


(32) 


Now, 


i  ( 


0.011 


HmA)  _ 
E(keV)3/2  r2 


(33) 


where  I(mA)  is  the  current  in  niA,  E(keV)  is  the  energy  in  keV,  and  r  is 
beam  radius  in  r^  units.  Thus, 


O'r) 


s2/3 


GC  -  r  Qin  t  I Ml! 

rq  1  k  E(keV)  T4/3 


and 


v  =  .152 


mj. 


1/3 


E(keV)1/2r?/3 


(34) 


(35) 


For  the  8.7  keV  case,  I(mA)  =  500.  The  background  plasma  is  characterized  by 

„  =  1.78  x  106  sec'1  when  n  =  105  cm'3. 

P  ® 

Since  '.r^  is  not  constant  for  the  cold  beam  case,  a  numerical 
calculation  of  beam  spreading  has  been  performed  using  the  Monte  Carlo  code 
described  in  Section  2.4.  The  results  are  shown  in  Figure  5  and*  compared  to 
the  warm  beam  results  of  Figures  1-3,  suggest  that  the  wave-particle  interaction 
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region  extends  much  deeper  into  the  atmosphere. 

While  the  computed  large  beam  spread  indicates  the  importance  of 
beai.i  plasma  turbulence,  the  results  are  not  completely  self-consistent.  The 
appreciable  beam  expansion  means  ■<  is  greatly  reduced,  so  that  it  is  un¬ 
likely  that  the  beam  can  remain  cold  at  lower  E-region  altitudes.  Thus 
Equations  (29)  and  (30)  should  be  replaced  with  warm  beam  expressions  some¬ 
where  above  100  kin.  The  intersection  of  (  and  v  occurs  at  low  altitudes,  but 
there  is  a  problem  with  performing  the  given  calculation  down  to  low  altitudes 
due  to  energy  loss.  Significant  loss  will  occur  below  120  km,  and  this  is  not 
accounted  for  in  the  calculation.  Thus  our  results  should  serve  simply  as  an 
indication  of  the  potential  of  turbulence  in  spreading  the  beam  under  discussion 


Sea  Li  ok  4 

Di£TA  1 L E U  ENERGY  DEPOSITION  AND  OPTICAL  EMISSION  PROPERTIES  USING  CODE  MCBE 

4.i  Computational  Model  a nd  P r ob 1 em  Parameters. 

In  the  previous  section,  we  presented  results  basically  analytical 
in  nature  giving  the  mean  radius  of  selected  e-beams  following  spreading  due  to 
turbulence  and  particle-particle  interactions.  The  contribution  from  the  latter 
effect  was  determined  by  using  the  analytical  result  of  Lin  and  Strickland  (1981) 
shown  as  the  vertical  line  in  Figure  4.  In  this  section,  we  augment  these  ana¬ 
lytical  results  with  more  detailed  results  from  code  MCBE.  These  will  in¬ 
clude  2D  (r,z)  contours  of  energy  deposition  and  various  other  related  quantities 
such  as  the  mean  radius  of  energy  deposition. 

Results  will  be  shown  for  particle-particle  interactions  alone  as 
we  1 1  as  for  the  combined  effects  of  wave-particle  and  particle-particle  inter¬ 
actions.  We  wish  to  make  it  clear  that  we  have  separated  the  regimes  for  these 
processes.  Thus,  MCBE  accounts  for  beam  spreading  by  wave-particle  interactions 
through  specification  of  the  radial  beam  profile  at  its  upper  altitude  boundary. 
The  beam  profile  is  calculated  by  the  model  specified  in  Section  2  which  is  im¬ 
plemented  in  the  code  WPDIF. 

Code  MCBE  was  developed  for  AFGL  under  a  previous  contract.  The  MCBE 
code  describes  the  transport  of  electrons  in  the  earth's  atmosphere  and  includes 
the  effects  of  the  earth's  magnetic  field.  The  code  solves  the  transport  problem 
using  an  extension  developed  by  D.  L.  Lin  of  the  Monte  Carlo  simulation  of 
multiple  scattering  originally  developed  by  Berger  (1963;. 
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Results  are  to  be  presented  for  the  two  beam  energies  8.7  keV 
and  38  keV.  The  initial  radial  distribution  for  the  8.7  keV 
beam  is  shown  on  the  top  graph  of  Figure  6.  Initial  distribution  1  is  the 
distribution  for  the  case  neglecting  wave  particle  effects.  Initial  distri¬ 
bution  2  is  the  distribution  for  the  case  including  the  effects  of  wave- 
particle  interactions.  Other  parameters  pertinent  to  the  8.7  keV  MCBE 
calculations  are  summarized  in  Table  3  below. 


Neglecting  Wave-  ! 

Particl  e  Effects  > 

Including  Wave- 
Particle  Effects 

1 

i Di stribution 

i 

1 

2 

!  Initial  A1 titude 

260  km 

120  km 

;  Initial  Angle 

70'  “rt  Searth 

70"  wrt  £ 

ea  rth 

'^earth1 

l 

. 5  gauss 

.  jL_ 

.5  gauss 

Table  3.  Summary  of  Input  Parameters  for  8.7  keV  Beam 


The  initial  radial  distribution  for  the  38  keV  beam  is 
shown  on  the  bottom  graph  of  Figure  6.  Initial  distribution  3  is  the  distri¬ 
bution  for  the  case  neglecting  wave-particle  effects.  Initial  distribution  4 
is  the  distribution  for  the  case  including  the  effects  of  wave-particl e  inter¬ 
actions.  Other  parameters  pertinent  to  the  38  keV  MCBjE  calculations  are  sum¬ 
marized  in  Table  4  below.  Since  the  38  keV  beam  is  assumed  to  be  warm  by  us. 
Equation  (16)  rather  than  WPOIF  was  used  to  generate  the  distribution  which  has 
been  placed  in  histogram  form  for  input  to  MC3E. 
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Neglecting  Wave- 
Particle  Effects 

Including  Wave- 
Particle  Effects 

Distribution 

Initial  Altitude 

Initial  Angle 

iB  ..  ! 

■  earth 

3 

212  km 

70°  wrt  3 

earth 

.5  gauss 

4 

165  km  j 

70°  wrt  §  +  ,  | 

earth 

.5  gauss 

Table  4.  Summary  of  the  Input  Parameters  for  the  38  keV  Bearn 


4.2  MCBE  Results. 

In  this  section  we  present  selected  results  from  the  MCBE 
calculations.  For  the  38  keV  beam,  the  assumption  that  wave-particle 
effects  and  particle-particle  effects  dominate  in  different  regimes 
should  be  valid.  The  region  in  which  wave-particle  effects  are  pre¬ 
dicted  to  dominate  is  above  160  km  and  the  energy  deposition  due  to 
particle-particle  effects  is  negligible  at  that  altitude. 

The  contour  plot  of  the  energy  deposition  D(z,r)[eV/cm  -e]  for 
the  38  keV  beam  is  shown  in  Figure  7.  For  the  case  which  includes  wave-particle 
interactions,  the  contour  line  for  a  given  level  extends  slightly  further  out 
in  radius  and  the  peak  occurs  (within  the  statistical  uncertainty  of  the  re¬ 
sults)  at  the  same  altitude.  Shown  in  Figure  8  is  a  plot  of  the  mean  radius 
of  energy  deposition  as  a  function  of  altitude.  The  calculation  which  includes 
wave-particle  effects  shows  an  overall  larger  mean  radius  with  a  maximum  dif¬ 
ference  of  less  than  20.  at  approximately  89  kilometers.  Figures  9  and  10 
show  the  column  energy  deposition  versus  tangent  point  distance  at  91.5  and 
115  km  for  trie  38  keV  beam.  The  Column  deposition  D  (?,r^  )  [eV/cm'’  -  ej  is 
obtained  by  an  integration  through  and  perpendicular  to  the  beam.  We  mav  thus 
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tude  (km) 


Energy  Deposition  (ev/cm  -e) 


Tangent  Point  Distance  (meters) 


Figure  10.  Column  Energy  Deposition  Versus  Tangent 
Point  Distance  at  a  Altitude  of  115  km 
for  a  38  keV  Beam. 
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view  this  quantity  as  being  projected  onto  a  plane  parallel  to  the  beam  axis. 
In  the  next  subsection,  we  will  discuss  how  the  column  deposition  may  be  re¬ 
lated  to  various  optical  emission  features.  Upon  transforming  deposition  to 
emission,  results  like  those  in  Figures  9  and  10  give  brightness  of  a  beam 
streak  in  two  dimensions. 


For  the  case  of  the  8.7  keV  beam,  there  is  an  overlap  between  the 
altitude  regimes  in  which  the  wave-particle  and  the  particle-particle  effects 
are  important.  The  WPDIF  calculation  for  this  case  indicates  that  the  altitude 
region  for  which  wave-particle  effects  are  important  extends  below  110  km. 

There  is  substantial  energy  deposition  for  particle-particle  effects  above 
110  km,  however,  and  it  was  necessary  to  seek  a  compromise  in  order  to  apply 
the  techniques  to  this  case.  An  altitude  of  120  km  was  chosen  as  the  altitude 
at  which  to  terminate  the  WPDIF  calculation  and  begin  the  MCBE  calculation. 

Shown  in  Figure  11  is  a  contour  plot  of  energy  deposition  for  the 

8.7  keV  beam.  From  this  plot,  one  can  see  that  the  energy  is  deposited  over 

a  much  larger  volume  for  the  case  in  which  wave-particle  interactions  are  included. 

_g 

Note  the  absence  of  the  10  contour  level  for  this  case. 

Figure  12  shows  that,  the  mean  radius  of  energy  deposition  is 
about  a  factor  of  four  greater  for  the  case  in  which  wa ve-pa rticl e  interactions 
are  included.  Figures  13  and  14  show  the  column  energy  deposition  versus 
tangent  point  distance  at  105  and  115  km  for  the  8.7  keV  beam. 


4.3  Op  t i ca 1  P rope r  t i es . 

In  this  sub-section,  we  relate  the  column  energy  deposition  rate 
such  as  shown  in  Figures  9,  10,  13,  and  14  to  optical  emission  properties. 
This  will  enable  us  to  estimate  the  brightness  of  e-beam  produced  optical 
streaks  as  might  be  seen  oy  a  distant  optical  imaging  system.  The  reference 
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Figure  11.  Contour  Plots  of  Energy  Deposition  for  a  8.7 

3 

keV  Beam  (units  are  eV/cm  -e) 
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Tangent  Point  Distance  (meters) 


Figure  13.  Column  Energy  Deposition  Versus  Tangent 
Point  Distance  at  an  Altitude  of  115  km 
for  an  8.7  keV  Beam. 


material  for  the  optical  aspects  of  the  discussion  to  follow  come  from 
Vallance  Jones  (1974)  and  O'Neil  et  al.  (1978). 


EMISSION  FEATURES 

We  have  chosen  to  examine  the  spectral  region  between  ■>.  3200  A 
and  8500  A.  Table  5  provides  a  list  of  some  of  the  more  prominent  features 
in  this  range  with  estimates  of  their  strengths  relative  to  that  of  N^+  IN 
3914  A  (Vallance  Jones  (1974)).  The  spectrum  becomes  more  complicated  at 
the  longer  wavelengths  which  is  the  reason  we  have  not  included  features  in 
the  table  for  such  wavelengths. 


EMISSION  FACTORS 


ship: 


We  designate  the  emission  factor  f  through  the  following  relation¬ 


al  (z  .r^  )  =  6.25  x  1012  fJD.(z,rJL  )  Rayleighs  (106  ph/cm2  -  s  -  4-  sr) 


where  4~I  is  the  column  emission  rate,  J  is  the  starting  beam  current  in  amps, 
and  the  constant  provides  the  proper  scaling  to  give  Rayleighs. 


Examples  of  f-values  are  given  in  Table  6.  Values  are  shown  for 
individual  prominent  features,  partial  and  entire  band  systems,  and  the  total 
of  the  systems  considered  for  wavelengths  less  than  8500  A.  The  given  value 
for  3914  A  provides  an  efficiency  of  5.4  x  10  (energy  emitted/energy  deposited). 
Assuming  10.  of  the  incident  beam  energy  is  backscattered ,  this  number  gives  an 
efficiency  relative  to  the  incident  beam  energy  of  4.9  x  10*^  which  is  close  to 
the  value  obtained  by  O'Neil,  et  al.  (1978). 


'(A) 

feature 

EMISSION  STRENGTH 

3370 

N2  2P  (0,0) 

30 

3536 

(1,2) 

1 

7  ! 

3576 

(0,1) 

20 

3582 

N2  IN  (1,0) 

7  ! 

1 

3709 

Np  2P  (2,4) 

2 

3754 

Np  IN  (1,3) 

6 

3804 

(0,2) 

8  | 

3884 

(1,1) 

4  ! 

1 

3914 

(0,0) 

100 

3997 

N2  2P  (1,4) 

4 

4058 

(0,3) 

;  3 

1 

4236 

N2  IN  (1,2) 

!  4 

4278 

(0,1) 

i  30 

! 

4709 

(0,2) 

6 

5577 

01  [lS  -  *D) 

100 

! 

6300 

01  {lD  -  3P) 

i  VARIABLE 

Table  5.  Prominent  UV  to  Middle  UV  Features 
And  Their  Emission  Strengths  Relative  To 
100  Units  Of  N*  IN  3914  A 


FEATURE 

EMISSION  FACTOR 

3914  A 

1.7  (-  3) 

N*  IN  SYSTEM 

2.6  (-  3) 

01  5577  A 

1.7  (-  3) 

2 P  3370  A 

5.6  (-  4) 

2P  3576  A 

3.3  (-  4) 

N?  2P  SYSTEM 

1.9  (-  3) 

n2  ip  system 

5.0  (-  3) 

(<8500  A) 

N2  vk  system 

8.3  (-  4) 

n*  m 

2.8  (-  3) 

(<8500  A) 

TOTAL 

1.4  (-  2) 

(SYSTEMS  plus 

01  5577  A) 

Table  6.  Approximate  Emission  Factors  For 
Prominent  Lines,  Bands,  And  Band  Systems 
Within  The  Interval  From  '\??0Q  A  To  8500  A 
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The  3914  A  f-value  just  discussed  was  obtained  from: 


1  0.8[N2]  1 

f 3  914  "  34  0\W2T^T0P“+  O'rSrol  14 


where  34  is  the  applied  eV/ion  pair  value,  the  weighted  density  ratio  gives 
the  relative  N^+  ionization,  and  14  is  the  applied  ratio  of  N2+  lons  t0  3914  A 
photons.  The  formula  fias  been  applied  to  the  100  km  region  where  [0]  is  unim¬ 
portant. 


The  emission  factors  other  than  for  3914  A  in  Table  6  were  obtained 
using  Tables  4.5  arid  4.9  -  4.13  from  Vallance  Jones  (1974).  All  such  values  are 
relative  to  f^g^. 

It  is  worthwhile  estimating  the  efficiency  of  all  features  considered 
in  Table  6.  To  do  so,  let  us  assume  an  average  wavelength  of  5000  A  which  cor¬ 
responds  to  an  energy  of  2.5  eV.  This  yields  an  efficiency  of  3.87  (energy/ 
incident  beam  energy)  allowing  from  10%  energy  backscatter.  This  can  be  com¬ 
pared  with  experiment  using  the  data  of  O'Neil,  et  al.  (1978).  Given  that  the 
beam  power  was  2  kW  we  estimate  an  efficiency  of  6%  from  their  Figure  9.  Con¬ 
sidering  all  uncertainties  and  the  fact  that  the  calculated  efficiency  was  de¬ 
termined  for  fewer  features  than  present  in  the  measurement,  the  agreement  in 
efficiencies  is  sa t i sfactory . 


ENERGY  DEPOSIT  ION  PROPERTIES 

To  utilize  tne  emission  factors  in  Table  6,  we  need  values  of  the 
i.ij'-"  cm  r-.-y  deposition  0  ( z. ,  r^  )  as  introduced  in  the  previous  sub-section. 

1  ;gur>>  if.  snows  tne  goom.-try  which  applies  to  the  determination  of  this  quantity. 
T  '■  1  ‘  ••i'o.i  mtooru '  i*.n  through  the  deposition  region  is 


40 


Dcl2-rJ 


)  * 


D(z,r(s))  ds  eV/cm  -e 


(38) 


where  d'.  is  an  element  of  path  length  along  the  1  ine-of-sight. 

EMISSION  KATES 

We  now  apply  the  information  in  Table  6  to  some  of  our  previous 
results  for  an  estimation  of  brightness  of  an  optical  streak  produced  by  an 
e-beam.  We  choose  to  consider  the  38  keV  beam  and  thus  address  ourselves 
to  Figures  9  and  10.  Attention  will  be  given  to  the  results  for  wave-particle 
and  particle-particle  effects  in  terms  of  the  sum  of  features  listed  in 
Table  6.  The  beam  current  J  in  Equation  (36)  is  taken  to  be  50  mA  since  this 
is  close  to  the  ECHO  IV  current  for  the  38  keV  shot. 


From  Equation  (36),  we  then  have 


4-rI  (z.r^ 


4.4  x  10*  Dc(z,r^  ) 


(39) 


Figure  16  shows  4 n I  as  a  function  of  r  at  several  altitudes.  Emission 
rates  such  as  these  can  be  used  directly  to  determine  the  requirements  of  optical 
imaging  systems  being  used  to  determine  streak  widths.  The  collection  rate  of 
Photons  by  the  system  is  given  by: 


F 


// 


4  I  ( z , ) dzdr^ 


photons/s 


(40  1 


A? 


where  is  the  collector  area,  R  is  the  detector-to-streak  distance,  and 
:.z  and  '.r^  give  the  surface  area  of  the  streak  viewed  by  the  detector.  As 
an  example,  consider  the  following  parameter  values: 


•  R  =  120  km 

•  z  =  Ar  =  10m 

•  Ap  =  100  cm^. 


For  the  brightnesses  given  in  Figure  9,  this  gives  for  F,  the  rather  weak 
value  of  31  photon/s  from  the  brightest  part  of  the  streak.  For  the  Davis 
observation,  this  value  should  increase  almost  an  order  of  magnitude  since 
the  beam  current  in  that  experiment  was  larger  by  a  corresponding  amount 
compared  to  that  of  the  ECHO  IV  38  keV  shot. 
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Appendix  A 


DERIVATION  OF  LINEAR  GROWTH  RATE 


The  best  developed  branch  of  plasma  physics  is 
linear  instability  theory  because  of  its  mathematical 
tractability .  Linear  instability  theory  can  identify 
those  modes  in  a  plasma  that  will  go  unstable  and  gives 
the  frequency  and  vector  wave  number  of  the  mode.  The 
theory  of  plasma  turbulence  is  the  attempt  to  understand 
how  these  linearly  unstable  modes  saturate.  In  this 
section,  we  develop  the  appropriate  linear  dispersion 
relation  for  this  problem  and  quote  some  of  the  pertinent 
results . 


The  pertinent  plasma  linear  instability  theory 
for  this  report  starts  with  the  Valsov  equation  and 
linearizes  about  a  perturbed  electric  field.  The  relevant 
equations  for  this  problem  are: 
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where : 


E 


B 


e  . 

J 

M  . 

J 

N  . 
J 


'Oj 


probability  distribution  function 
for  jth  species  as  a  function  of 
position,  velocity,  and  time; 

electric  field; 


magnetic  field; 

t  h 

charge  of  the  j  species; 

,  , .  .  th 

mass  of  the  j  species; 

j1*1  species  particle  density  =  J~ dj^f  ^  ; 

unperturbed  j**1  species  particle 
density;  and 

species  -  neutral  collision 
frequency . 


The  righthand  side  of  Equation  (A.l)  is  the  so-called  BGK 
collision  term  (Parley  (1963))  and  has  been  demonstrated  to 
adequately  represent  the  effect  of  charge  particle-neutral 
collisions  on  linear  growth  rates.  Of  course,  Equation  (A. 2) 
is  just  Poisson's  equation. 

Equations  (A.l  )  and  (A. 2)  are  linearized  by 

assuming : 


B(_x,t) 

=  B0(x,t)z; 

E(_x,  t  ) 

~  i  A 

=  kE(k,-  )e  ‘  (electrostatic  approximation) 

f j( ) 

-  fj0<»)  + 

.  iA 
'Oe  ; 

,  ,2 

>1 

o 

l< 

i 

l< 

o 

c_j. 
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2kBTi  2 

— —  =  2v^.  =  thermal  velocity  squared 

j 

( a)  -  k  •  *  )  ; 

radial  wave  frequency; 
wave  number; 

"t  h 

drift  velocity  of  j  species. 

To  make  this  standard  electrostatic  perturbation  fit  the 
problem  at  hand,  we  also  assume  we  have  three  species: 

(1)  plasma  ions  (i), 

(2)  plasma  electrons  (e),  and 

(3)  beam  electrons  (B). 


The  beam  electrons  have  a  net  motion  in  the  Z  direction: 


V  = 

-Oi 

0; 

v  - 

-Oe 

0; 

V 

V 

-OB 

Z 
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Under  these  assumptions,  the  first  order  equation  from 
Equation  ( A  -  1 )  can  be  written: 


D  „  ,  ,  iA 

iK  f-(k.v.u)e 


_  "  [tf1  E».')'0jO-)  -  "j  fd  -iA 


where : 


D 

Dt 


at  +  -v  *  Tx  +  m1  (-  x  5)  Tv  +  * 

3 


J 


(  A.  A) 


Equation  (A. 2)  can  be  solved  for  by  integrating 

over  the  unperturbed  orbit  defined  by  the  operator, 

Equation  (A.l).  f.(k,v,u>)  can  then  be  integrated  over 

J 

velocity  as  prescribed  by  Equation  (A. 2)  to  obtain  a  linear 
dispersion  relation.  This  procedure  is  well  documented  in 
several  sources,  for  example,  Clemmow  and  Dougherty  (1969)  and 
St  ix  (19(52).  The  unal  result  ol'  this  procedure  is: 


1  +  Z  H .  (u,k)  =  o 

J  J 


(A.  5) 


Hj(ft).k) 


W_Bi_  ^  +  iWjGj(u),k) 


,2 02  1  -  v.G  (u>,k) 

k  vthj  J  J 


(  A  .  6 ) 


Gj  (  u),  k  ) 


-  i 


"^vthj  kZ 


-  X  . 
,  J 


y  I  (  X  .  )Z(  ;  .  )  ;  (A  ■ 
Z,  nv  j'  jny  ’  *'• 


n  =  -oo 


r>n 


where : 


2 

0) 

PJ 

= 

4irn.e./M.  =  plasma  frequency; 

J  v)  J 

♦ 

= 

w  +  iVj  -  kzvz; 

;nj 

= 

l/y  2  _  nftj)  thermal  velocity 

Xj 

= 

(,thj  /"j  ki>2; 

kx 

= 

k  vector  perpendicular  to  B-field; 

kz 

= 

k  vector  parallel  to  B-field; 

Xn 

= 

t  h 

nx  modified  Bessel  function;  and 

Z 

s 

plasma  dispersion  function 

Fried  and  Conte  (I960). 


G.  [Equation  (a.  7)3  is  known  as  the  Gordeyev 
J 

integral ,  Cl emmow  and  Dou^hertv  (1969).  To  extract  a  mean  i  iu»  t'u  1 
amount  ol  information  from  Equation  (A. 5),  the  usual  "cold" 
miasma  limit  is  taken.  This  limit  is: 


\  .  <<  1 
J 


\jn 


>>  1 


u ) 


PJ 


>>  V  . 

J 


(  A.H) 
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which  is  equivalent  to  saying  kv^j  <<  uo  ,  ,  and 

collisions  are  negligible.  In  this  limit,  only  the 
n  =  -1,0, +1  terms  in  Equation  (A. 6  )  are  retained  and 
the  Bessel  function  and  plasma  dispersion  function  are 
expanded  by  the  following  prescription  (  Abramowitz  and 
SI  I'l'iin  (1  !>()■!  )). 


I  (X) 
n 


1  ' 

(2X)n 


Z(?) 


These  expansions  led 
j  =  i  ,  e  : 


to  the  following  result  for  H.  when 

0 


H.(u,k) 


(  A.O 


The  limit  [Equation  (A. 8,1  is,  however,  not  valid  for  the 
beam  electrons.  In  this  case,  it  will  turn  out,  the 
electrostatic  perturbations  of  interest  will  process  the 
f o 1  I ow i n g  proper! y : 


'Bu 


1  VB  -  kZVZ 


M  , 

/J  J 


1 


(  A  .  1 1  ) 


h? 


BiBWaaii 


w 


For  this  scaling,  the  Landau  prescription  for 
expanding  the  plasma  dispersion  function  is  operative 
(Slix  (.1964)),  and  leads  to  the  result: 


HgU.k)  = 


ie 


B 


=  -i 


'Z 

3 


U). 


PB 


f0B(^/kz) 

nB 


(  A  .  1  1  ) 


Combining  the  approximations  of  Equations  (A. 9)  and  (A.ll) 

into  the  plasma  dispersion  relation,  Equation  (A.E),  and 

2  2 

neglecting  the  ion  contribution  ( to  >>  w  we  obtain: 

&  &  pe  pj 7 ’ 


1 


i  c 


B 


(A.  12) 


Equation  (A. 12)  is  the  standard  approximate  dispersion 
relation  for  this  two  stream  or  bump  on  tail  instability. 
When  this  approximation  is  legitimate  and  the  growth  rate 
is  much  less  than  the  real  part  of  the  frequency: 

w  =  wR  +  iy,  >>  y 

where  y  =  growth  rate,  Equation  (  A.  12)  can  readily  be 
solved  bv  first  setting  y.e^  to  zero  and  solving  for  u)n ; 
then  expanding  the  right-hand  side  in  terms  of  \  and 
equating  it  to  Eg. 

The  result  of  these  manipulations  are: 
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QU) 


(A.  14) 


.  2  n2  w  2  2  2. 
w(u)  - 0  )(w  -S2 

2  r2  2 7“ 
( 2  ui  - fi  -  w  ) 


K2 


,  2  r.2w  2  2, 

(a  -•<  )(“  -<*>p) 

2  2  ,  2  2~ 
CJ  \  Cl)  *“  it  * 


(  A  .  1 ) 


where  we  have  to  drop  the  "e"  subscript  and  redefine 

a,  _  r  .  As  always,  it  is  argued  that  the  dominance  to 
K 

oscillation  will  be  the  one  with  the  largest  growth  rate. 
The  procedure  is  to  maximize  Q(oj)  given  an  .T  and  w  . 

i*' 

This  is  graphically  worked  out  in  Figure  17.  Assuming 
5  -3 

np  =  10  cm  and  =  0.58  gauss, we  have: 


U  = 


eB 

M  c 
e 


,n7  -1 

10  sec 


4-rrn  e 


- —  =  1.66  x  107  sec  1 


w  2 


=  2.75 


~  =  0.6  (see  Figure  17) 


=  6 . 0  x  10®  sec  ^  . 
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FIGURE  17 


hi  as  a  Function  of  w  /.Q  to  Maximize 

P 

Q(uO  in  Equation  A. 13 


(St ix  (1964)) 


Using  this  value  of  ui  in  Equation  (A. 13) 

gives : 

Q  ,  (ui)  =  0.4 

^ma  x  ■  ' 

Next,  is  maximized  by  choosing  an  appropriate  k„ 

D  L 

[see  Equation  (A. 13  )] .  This  is  approximately  accomplished 
by  setting: 


v 


7 


thR 


(  A.  16) 


In  the  limit,  <<  v^,  Equation  (A. 1(5)  can  be  used  in 

Equation  (A. 11)  to  transform  Equation  (A. 13)  into: 


Y  =  0.76  Q 


(  A  .  1  7  ) 


where  a  =  nD/n  .  Of  course,  the  ratio  k,/k„  is  determined 

d  C  i-  L 

by  Equation  (A.  13)  once  is  chosen  to  maximize  Q(u.').  This 

ratio  is  shown  in  Figure  18  as  a  function  of  (u  / E )  •  For 

P 

our  analysis: 


which  translates  into  an  angle  of  Gj  =  45°,  where  0j  is 
the  angle  between  k  and  B. 


The  question  this  analysis  does  not  answer  is 
the  effect  of  collisions  on  the  instability.  Roughly 
speaking,  we  can  see  the  effect  in  the  expansion  condition 
Equation  (  A. 10)-  With  negligible  collisions,  the  expansion 
(A .  11  )  was  with  respect  to: 


UR  +  iY  +  kZVZ 
kZ  VthB 


<<  1 


Obviously,  if  we  substituted  iv„  =  iy  in  the  computation 

r> 

lor  Hg(o. ,k),  we  would  get  the  same  numerical  answer. 

Thus,  a  quick  and  dirty  approximation  such  that  collision 
frequencies  cut  off  this  instability  would  be  the  condition 

Y  =  Vn  •  (A. 18  ) 
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<N  H  <N  N  0.06 
it  it 


0.006  - 
0.004  - 


0.002 


(u  /n>‘ 


FIGURE  IE.  Ratio  of  Perpendicular  to  Parallel 
Electric  Field  that  Corresponds  to  Maximum 
Q  in  Equation  (A. 14) 

(Slix  (1964)) 


This  is  really  an  argument  that  the  instability  will 
stabilize  with  the  growth  rate  which  is  of  the  same  order 
as  the  collision  frequency.  This  completely  neglects  the 
contribution  of  the  denominators  in  Equation  (A . 7  )  .  A 
more  precise  statement  would  require  the  numerical  solution 
of  the  exact  dispersion  relation,  which  is  a  reasonably 
tractable  computational  problem. 


Appendix  B 


DEFINITION  OF  COHERENCE  TIME 


To  explain  the  diffusion  of  auroral  streaks, 
we  hypothesize  the  linear  unstable  waves  introduced  in 
Section  3  grow  and  saturate  to  some  level  of  turbulence, 
producing  an  E-field  that  makes  the  particles  E  x  B 
drift  perpendicular  to  the  axis  of  the  beam.  The 
turbulence  is  pictured  as  separate  regions  of  coherence, 
where  within  each  region  a  definite  electrostatic  wave 
eixsts  with  a  particular  w  and  k  nominally  of  the  value 
predicted  by  the  linear  theory.  If  the  beam  particles 
traverse  these  regions  experiencing  a  reasonable  constant 
E-field  per  region,  the  guiding  center  of  the  particles 
will  be  displaced: 


Ax  =  v 


drift  Tcob 


where : 


drift 


E  x  B 


cB 


x 

coh 


time  particle  transverses 
coherence  region;  and 


c 


velocity  of  light. 


In  terms  of  gyro-radius,  this  Ax  can  be  expressed  as: 


rg  MvBsinOB  T  coh  '  (H,l) 

Using  this  random  displacement,  the  radial  diffusion  of 
the  electron  beam  can  readily  be  computed.  The  objective 
of  turbulence  theory  is  to  calculate  the  coherence  time 
and  the  magnitude  of  the  E  -field. 

COHERENCE  TIME 


To  estimate  coherence  time,  we  estimate  the 
length  of  the  coherence  region.  To  see  how  this  works, 
consider  Figure  IS.  Part  (a)  of  Figure  19  is  a  picture  of 
a  spatially  uniform  wave.  From  Section  2,  we  know  this 
wave  has  wave  fronts  that  move  with  velocity (w/k )  =  v^, 
the  phase  velocity.  Part  (b)  shows  a  wave  packet  con¬ 
taining  the  wave  in  Part  (a).  The  wave  fronts  in  the  wave 
packet  still  travel  at  the  phase  velocity,  but  the  wave 
packet  itself  will  travel  of  the  group  velocity.  A  beam 
particle  transversing  the  packet  at  vz  will  be  in  step 
with  a  wave  front,  and  the  time  within  the  packet  will  be: 


T 

coh 


L  =  coherence  region  length 
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(a) 


FIGURE  IS (a).  Plane  Wave  with  Phase  Velocity  v  ^ 


I 

FIGURE  19(b).  Wave  Packet  with  Phase  Velocity  vz 
and  Group  Velocity 

In  both  pictures,  electron  velocity 
is  phase  velocity;  therefore,  electrons 
see  constant  E-field. 


If  this  picture  is  valid,  the  wave  packet  will  have  to  be 
at  least  one  wave  front  wide.  Without  further  justifica¬ 
tion,  we  estimate: 

Tcoh  =  A3  X  |vi"  I  ’  A3  *  1  III. 2) 

i  Z  g  I 

where  is  an  arbitrary  constant  introduced  as  a  measure 
of  the  uncertainty  in  this  estimate. 
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From  Equation  (3.12),  we  can  show  the  group  velocity 


along  the  B-field 
so  rc0^  will  be: 

will  be  about  0.5  (O'Neil 

and  Walmber 

Tcoh  = 

A  — 

3  w 

(l).:i  ) 

Note  that  this  argument  justifies  the  statement 

the  beam  particles  will  see  a  constant  E-field 
a  coherence  region  for  any  group  velocity. 

that 

within 

E-FIELD 

SATURATION 

To  estimate  the  E-field  saturation,  we  use  two 

assumpt ions : 

Ttr 

A4Y_1  ,  A4  <  1; 

(  B  .  1  ) 

E2/8tt  = 

A5nB^B^B  ’  A5  ~  1, 

(B.5) 

where : 

-2 

Ttr 

kzeE^/M,  and 

A4’A5 

arbitrary  constants. 

In  Equation  (  13.4),  tt  is  referred  to  as  the  trapping 
-1  1  r 

time.  t.  =  u0  is  the  bounce  frequency  of  electrons 

tl  jD 

trapped  near  the  bottom  of  a  spatially  oscillating  electric 
f ield . 
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A  linearly  unstable  wave  grows  in  a  plasma 
because  a  reasonable  number  of  beam  particles  undergo  an 
osoil juiion  in  phase  with  each  other.  When  the  trapping 
time  of  these  particles  gets  to  be  of  the  order  of  the 
slowest  timescale  of  the  linear  instability  analysis 
(the  growth  rate),  the  zero  order  orbits  CEquation  (B.l)j 
are  modified  to  the  extent  the  analysis  breaks  down, 
hence.  Equation  (B.  1  ). 


Equation  (B.h)  is  a  statement  that  the  energy 
associated  with  the  electric  field  will  scale  with  the 
heating  of  the  plasma  beam.  The  more  wave  fluctuations, 
the  more  heating  of  the  plasma  beam. 

The  assumptions,  Equations  (if -  1 )  and  (B..1),  are 
difficult  to  rigorously  justify.  They  do  reproduce  the 
scaling  of  several  1-D  computer  simulations  and  are  not 
contradicted  by  weak  turbulence  theory  (Davis  ot  al.  (1971) 
and  Ha  li  inn  i)  ot  al.  (  1978)).  The  advantage  of  these  assumptions 
i  that  eombiiU'd  with  Equation  (A.  17),  they  form  a  closed  set 
! or  t He  rnlculat  ion  of  the  saturated  E-field. 

Starting  with  Equation  (  C.l),  we  have: 


Letting  k7  -  u./vz,  we  can  rewrite  Equation  (  B.(5)  to 
read  : 


IT 

8  r" 


,  2 


2) _ 

Z  /  2  . 

/  cos  Gj 


f>? 


4 


1_ 

(x 


Then,  using  Equation  (6.5),  we  can  rewrite  the  above  to 
read : 


nB  2 

X  tfZ 

nBkBT 


2  2 

w  co  u  0 
e  2 

- 4—  COS 

Y 


(IS.  7  ) 


Equations  (6.1)  and  (B.7)  can  now  be  substituted  in 
Equation  (A.  17)  to  obtain: 


(y/a4  )‘ 


=  T 


-  ? 
tr 


[- 


n c.  Oa  2  2  2  2 

76  12Q  a  u)  w  cos 
e 


2/5 

1  A/i  J 


(  13.  F) 


This  gives  the  bounce  frequency  as  a  function  of  parameters 
derived  from  linear  plasma  theory. 


BEAM  DIFFUSION 


Equation  (13.8)  can  now  be  used  to  calculate  the 

random  step  size  of  the  guiding  center  of  beam  electrons. 

From  the  definition  of  ,  we  have 

tr 


(  13.6) 


Substituting  the  proper  angle  dependence  into  Equation  (:3.9) 
and  using  Equations  ([;.]),  (B.3),  and  (B.8).  we  can  derive: 


,  tan  Or  • 

A_x  _  . _ _1  4  tt 

r  3  tan  Qn  2  2 

g  B  <*>  tTr 
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A_* 

r 


rA5 

A3l4 


4  it  tan  L 


i2/5 


tan  & 


-  4j  0.7G  HQ  «2a;2ajg  cos2  G  j  .  ( 

B  u'  *-  J 


Assuming  the  following  parameter  values  with  the  Echo  IV  experiment  in 
mind, 


'I 


Q  = 


(  n  .  l  i) 


we  have: 


L  x 

r 


17.6  A 


2/5 


4/5 


3  V  A 


Then,  using  Equation  (B.I),  we  can  reduce  the  random  step 
size  expression  to: 


17.6 


where  I  is  in  amps  and  is  in  keV. 


(13.  12) 


ll .  10) 
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COMPARISON  WITH  EXPERIMENT 


Equations  (A.  17),  (B.K),  and  (B.12)  can  now  be 
used  to  compare  with  the  ECHO  IV  experiment.  To  make 
this  comparison,  we  first  set  our  arbitrary  constants 
A  A.,  and  Ac  equal  to  unitv.  These  constants  indicated 

.5  b 

the  sensitivity  of  the  estimates  we  made  for  coherence 
time,  trapping  time,  and  beam  thermal  spread,  respectively. 


First,  we  test  the  assumption  that  the  instability 
turns  off  at  150  kin  due  to  electron-neutral  collisions. 

Using  Equations  (B.«)  and  (B-H),  the  growth  rate  will 
equal : 


y 


(6.8  x  106) 


-1 

sec 


.  (!’>.!  3 ) 


We  now  insert  into  rg  and  <q,  their  values  as  estimated  for  pulse  10  which 
gave  the  most  diffuse  optical  streak.  These  values  are  taken  to  be  r^  -  10 
gyro-radii  and  .g  =  .00032  which  yield: 

Y  =  4.3  x  104  sec  1 

According  to  Berger,  et  al .  (1970)  this  corresponds  to  a  collision 

frequency  at  an  altitude  of  120  km.  At  150  km,  the  collision  frequency 
3-1 

is  1.2  x  10  sec  ,  a  factor  of  36  lower. 
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The  second  assumption  we  would  like  to  test  is 
the  diffusion  attributed  to  the  random  step  size, 

Equation  (  [5 .IE).  Again,  consider  pulse  10.  1  or  tquation  (13.12)  with 

A 3  '■'a  Al,  i ,  we  have: 

/.  X  in  a.  4  / 5 

—  =  17.6  a 

r 

g 


Obviously,  as  the  beam  descends  through  the  atmosphere, 
ct  decreases  and  this  step  size  will  decrease.  An  over¬ 
estimate  of  the  beam  spreading  would  be  to  use  a  constant 
step  size  equal  to  the  initial  or  maximum  step  size.  In 
this  case,  the  maximum  beam  spread  will  satisfy  the 
inequality : 

ilB  <  =  /F (17.6  aj/5)  ,  (  li.3-0 

rg  rg 


where  N  is  the  number  of  random  steps.  This  can  be 
estimated  by: 


N 


fl ight  time 
~  coh 


1.2  x  103 


Thus,  an  upper  bound  for  streak  radius  is: 


r  +  r 
8 _ 


B 


r 

f- 


1.97 


which  if-  about  five  times  smaller  than  the  observed 


06 


streak  width. 


